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P-h ■ Abstract 

I> 

In other work, the generalized counting quantum Turing machine (GCQTM) was 

studied. For any N this machine enumerates the first 2 integers in succession as bi- 
nary strings. The generalization consists of associating a potential with read 1 steps 
only. The Landauer Resistance (LR) and band spectra were determined for the tight 
binding Hamiltonians associated with the GCQTM for energies below the potential 
height. Here these calculations are extended to energies both above and below the 
barrier height. For parameters and potentials in the electron region, The LR fluc- 
tuates rapidly between very high and very low values as a function of momentum. 

i-£h I The rapidity and extent of the fluctuations increases rapidly with increasing N. For 

N = 18, the largest value considered, the LR shows good transmission probability 
as a function of momentum with numerous holes of very high LR values present. 
This is true for energies both above and below the potential height. It is suggested 
that the main features of the LR can be explained by coherent superposition of the 
component waves reflected from or transmitted through or across the 2 Ar_1 poten- 
tials present in the distribution. If this explanation is correct, it provides a dramatic 

K> | illustration of the effects of quantum nonlocality. 



1 Introduction 



Since the early 1980s quantum computation has been under active develop- 
ment in several directions. This includes the development of efficient algo- 
rithms [1-3], simulation of physical systems[4,5], quantum error correction [6], 
and use of quantum gates [7]. In this work and in earlier work on Quantum 
Turing Machines [8,9], computations were assumed to proceed smoothly with 
no potentials present to cause reflections or decay of the transmitted state. 

The importance of including these potentials and the resulting degradation 
of performance has been emphasized by Landauer [10]. He has noted that 
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potentials can be present because of environmental influences and physical 
limitations and errors [11,12] which can occur in realistic models. 

In previous work [13,14] on quantum Turing machines (QTMs) potentials were 
introduced by associating them with different types of steps. Isolation from the 
environment is assumed. In this case the Feynman [15] Hamiltonian for these 
generalized quantum Turing machines (GQTM)s, can be also be described as 
a sum of tight binding Hamiltonians each of which has a different potential 
distribution along a one dimensional path of computation states. Which of 
these Hamiltonians applies to the GQTM is determined by the initial state. 

The effect of these potentials was studied in detail for the counting general- 
ized QTM (GQTM), which enumerates the nonnegative integers in succession 
as binary strings. The potential distribution associated with this GQTM, is 
closely related to the heirarchical sequence [13,16,17]. This sequence gives the 
number of Is appearing before the first in the successive iteration of integers 
represented as binary strings. 

Both the band spectra and the Landauer Resistance [18,19] were calculated as 
a function of energy for computations which enumerate the first 2 N integers 
as binary strings. The calculations were limited to energies below the height 
of the potential (bound state energies). 

In this paper the calculations made for the band spectra and LR for the 
counting GQTM will be extended by removal of the limitation of the energy 
range to energies below the barrier. Since the physical model used corresponds 
to head motion on a 1-D lattice of qubits with a Hamiltonian constructed from 
step operators, the full energy range remains finite. 

In the next section the physical model and main points of the description 
of the Hamiltonians for GQTMs will be reviewed. The use of the transfer 
matrix to calculate both the band spectra and LR will be briefly reviewed in 
Section 3. Following this the counting GQTM will be briefly described. Since 
the material has already been presented elsewhere, the discussion will be quite 
brief with additional material incorporated by reference. 

Results of calculations of the band spectra and the LR will be presented 
and discussed in Section 5. The LR, plotted as the log LR as a function of 
momentum, will be seen to fluctuate rapidly between very large and very 
small values. The fluctuations get more rapid and extreme as iV is increased 
and the LR increases as the potential height is increased. The association 
between the LR and band spectra for the associated periodic system obtained 
here is similar to that obtained by others [20,21] in that peaks are associated 
with gaps 

In Section 6 problems with extending the calculations to moderate values of 



iV are briefly discussed. The possible importance of coherent superposition 
of wave components reflected from the many potentials in the distribution in 
producing the observed fluctuations in the LR will be discussed. It is noted 
that this provides a good example of quantum nonlocality. 



2 Hamiltonians for GQTMs 

2. 1 The Physical Model 



The physical model for QTMs consists of a two-way infinite lattice of systems 
each of which can assume any state in a finite dimensional Hilbert space 7i s . 
If 7i s is two dimensional the systems are called qubits, However this term 
is used here irrespective of the dimensionality of Ti s . A head, described by 
a finite basis of states |Z) with leL, moves along the lattice interacting with 
the lattice systems at or adjacent to the head position. Possible elementary 
actions include any combination of head motion one site in either direction, 
lattice system state change, and head state change. 

The system states are all assumed to lie in a separable Hilbert space 7i. A par- 
ticular basis, the computation basis, defined by the set of states {|/, j, S)} and 
which spans 7i, is used. Here l,j refer to the internal state and lattice position 
of the head. The qubit lattice computation basis state \S) = <S>m=-oo\S( m )) 
where S is a function from the integers to the qubit state labels (e.g. {0, 1}) 
such that S(m) ^ for at most a finite number of values of m. This condition, 
the tail state condition, is one of many that can be imposed to keep the basis 
denumerable. 



2.2 Tight Binding Hamiltonians for GQTMs 



A bounded linear operator T is associated with each GQTM. T or its adjoint 
is referred to as a step operator for the GQTM because iteration of T (or its 
adjoint) corresponds to the successive steps in the forward (or backward) time 
direction of the GQTM. An infinitesimal time interval is associated with the 
steps of T so that it can be used directly to construct a time independent 
Hamiltonian according to Feynman's prescription [15]: 



H = K(2 - T - T f ) 



where K is a constant. Note that if T — Y where Y is the bilateral shift along 
the lattice, then H is the kinetic energy of free head motion on the lattice. 



As such it is equivalent to the symmetrized discrete version of the second 
derivative, (—h 2 /2m)d 2 /dx 2 . 

For each GQTM the associated step operator T is defined as a finite sum over 
elementary step operators. That is T = J2i s Tis where the sum is over all head 
state labels leL and all qubit state labels seS where S is the set of qubit state 
labels. For binary qubits S = {/, oo}. T} jS corresponds to the action taken by 
the GQTM associated with T when the head in state \l) sees or reads a qubit 
in state \s). 

For each I, s 7] jS = 7/ jS W/ jS where 7/ jS is a positive real constant < 1. Wi >3 is a 
sum over all lattice positions of products of projection operators, bit and head 
state change, and head position change operators. The specific definition of 
Wi yS appears elsewhere [13]. 

The generalization from QTMs to GQTMs consists in allowing values of 
7/ jS 7^ 1. Thus every QTM is a GQTM with 7; s = 1 for all l,s. Note that 
the generalization does not expand or change the definition of quantum corn- 
put ability. 

The main condition imposed on T and its adjoint is that they be distinct path 
generating in some basis B which spans H. This means that iterations of T or 
Tt on any basis state in B, generate a path of states that, up to normalization, 
are also states in B. The paths are distinct in that no two paths join, intersect, 
or branch. 

Although B can be any basis which spans 7i, (see [22] for details), unless 
explicitly specified otherwise, B is limited here to be the computation basis, 
{\l,j,S)}. A path in the B is defined as a set of states in B ordered by 
iteration of T or T*. If the state \l m ,j m ,S m ) = \m,i) is the mth state in 
some path i, then the m + 1st and m — lth states in the path are given 
respectively by \m + l,i) = W lmiSm (j m )\l m ,j m , S m ) = \l m +i, j m +i, S m+1 ) and 
\m-l,i) = ^ m _ uSm _ i(jm _^\l m , j m , S m ) = \l m -idm-i,S m -i). These equations 
must be modified if either T or T^ annihilate \m, i). 

If T is distinct path generating in B then T can be decomposed into a direct 
sum T = J^iTiPi where Pi projects out a subspace Hi of Ti spanned by a 
subset Bi of B. For each % iteration of T and T' on the states in Bi orders 
them into a path, denote the computation path i. The action of Tj and T- on 
each state \m,i) = \l m ,j m , S m (j m )) which is not a terminal path state is given 
by 



± |'mj Jm; bm) ~ ll m ,S m (j m ) \<-m+l, Jm+l, bm+l) 
-L Km) Jmi ^> m ) lflm-i,jm-i,Sm-i(jm-i)\'"rn— 1) Jm— li &m— 1/ ■ \*>) 



Additional details on the decomposition and its description as a sum of weighted 
shifts is given elsewhere [14]. 

Corresponding to the decomposition of T is a decomposition of the Hamilto- 
nian, Eq. 1, H — Ya HiPi where for each i, 

H l = K{2-T l -T}). (3) 

As is noted elsewhere [13] Hi can also be written in the form 

Hi = K{2 -Ui- U}) + Vi (4) 

where Ui is a shift along the path i. Ui and U\ are defined by Eqs. 2 with 
7 M = 1 for all I, s. 

K(2 — Ui — Uj) is, formally, the kinetic energy associated with the evolution 
of the computation along path i. For the physical model used here it is also 
the kinetic energy of head motion on the qubit lattice. (The kinetic energy 
depends on the magnitude but not the direction of the head momentum on 
the lattice.) 

The potential Vi = K{Ui — Tj + U\ — Tj) is a nearest neighbor off-diagonal 
potential with matrix elements on internal path states given by 

(m',i\Vi\m,i) = K[(l - 7i m ,s m (j m ))<*m',m+i + (1 - 7i m -i,s m _iCj m _i)<*m',m-i]-(5) 

In Schrodinger evolution under the action of H, Eq. 1, the choice of which Hi 
is active is determined by the initial state \l/(0). If ^(O) is a superposition of of 
computation basis states in one path i, then Hi is the only active component. If 
^(O) is a superposition of computation basis states in different paths, then all 
Hi are active for just those paths i in which \l/(0) has a nonzero component. In 
this case the states in the different paths evolve coherently which corresponds 
to quantum parallel computation [8]. 



3 Transfer Matrix 



For each computation path % one is interested in the evolution of the system 
state along the path under the action of Hi. For paths in which the region of 
nonzero potentials is finite in extent with potential free regions on both sides, 
the probability of completion of the computation is given by the transmission 
probability of the state through the potential region. This can be obtained by 



use of the transfer matrix following methods developed by Erdos and Herndon 

[19]. 

To be specific, assume the potential is nonzero in a region extending from 
state \a,i) to \a + m, i) and is zero elsewhere on an infinite path i. The states 
on either side of the region can be written in the form, 



a-l 



*!= Y, (^e ik J + Be- ik J)|j) 



j=-oo 

CO 



j=a+m+l 



(6) 



Here j is the path state label (the path label % in \j, i) is suppressed) and k is 
the momentum in the potential free regions. 

The form of \l/jj, the state in the region II of nonzero potential, depends on 
the form of the potential. In general \I/jj = X^=™(Cje llj -J + Dje _llj J)|j) where the 
momenta lj and coefficients Cj, Dj depend on the potential height, Eq, 5, at 
| j). Both the momenta and coefficients are constant over subregions of constant 
potential height. For any subregion for which the energy corresponding to the 
momentum lj is less than the potential height for the subregion, lj is replaced 
by ilj. Additional details are given elsewhere [19,14]. 

The four complex coefficients A,B,F,G are completely determined by bound- 
ary conditions for the case under study and the properties of the transfer ma- 
trix. For computations it is convenient to let Regions I and III be the initial 
and final state regions. In this case, the A, B, F coefficient terms of Eq. 6 are 
respectively the initial incoming state, the reflected term and the transmitted 
term. In this case G = as there is no incoming term in the final state region. 
Different boudary conditions are used if \& represents an eigenstate of Hi. 

The transfer matrix Z relates the coefficients F, G to A, B according to [19] 



(7) 



Z is a unimodular (determinant=l) 2x2 matrix which also satisfies [19] 
Zu = Z 22 

Z\2 = Z 21 . 



The above definitions can be iterated to define Z as a product of transfer 



matrices for individual subregions of constant potential. If there are M such 
regions then Z is the product of M 2 x 2 matrices, taken in the order in which 
the subregions occur in region II [19]. 

The Landauer Resistance (LR) [18,19], which measures the resistance to trans- 
mission through the potential region, is defined by LR = \B\ 2 /\F\ 2 with G = 0. 
For a computation on a path i, the probability of completion is given by 
1/(1 + LR) = \F\ 2 . The properties of Z can be used to obtain [19] 

LR=\Z\ 12 \ 2 (9) 

It is clear from the above that the LR depends both on the path and on the 
system energy (i.e. the energy of the head as it moves on the lattice). 

The transfer matrix is also used to determine for each path % the energy band 
spectra for the infinite 1-D crystal whose unit cell is given by the potential 
distribution of Hi, Eq. 5 (the associated periodic system). As is well known 
[23,20,21] the band and gap distribution is obtained from a plot of \TrZ\/2 as 
a function of momentum or energy. Regions in which \TrZ\/2 < 1 correspond 
to bands; regions with \TrZ\/2 > 1 correspond to gaps. 



4 The Counting GQTM 



An interesting example which can be used to illustrate the foregoing is the 
counting QTM [13,14]. Starting with a string of 0s this machine generates 
nonnegative integers as binary strings in succession by repetitve addition of 
1 mod 2. One marker locates the start of the strings on the qubit lattice and 
another marker, if present, locates the maximum length of the strings. 

A step operator T which carries this out can be given as a sum of 7 terms. 
A simple but nontrivial representation as a generalized QTM is obtained by 
setting, for all /, 7^ = 1 for all s ^ 1 and 7^ = 7 < 1. The resulting step 
operator is given by 



T= J2 (QoPojuPj + wQoPijuPj + QiPojuPj 

j=-oo 
+ wQ 1 P 2 jU^P j + ^Q 2 V xj Pijv)Pj + w ] Q 2 v xj Pq j uP j 

+ wQ 2 P 2j uP 3 ) (10) 

The projection operators Qi, P s ,j Pj refer to the head in state |Z), the site j 
qubit in state \s), and the head at site j. w is a shift mod 3 on the three head 
states (wQ m = Q m+ iw mod 3) and u shifts the head along the lattice by one 



site (uPj = Pj + \u). The need for markers is accounted for here by choosing 
the qubits in the lattice to be ternary with states |0), |1), |2). |2) is used as 
a marker and |0), |1) are used for binary strings. The qubit transformation 
operator v X j = a X j(P j + Py) + P 2 j exchanges the states |0), |1) and does 
nothing to the state |2) for the site j qubit. 

The adjoint T' is defined from T in the usual way noting that the operators 
for the head states, head position states, and qubit states commute with one 
another. Note that v X j = v X j. The explicit form of T* as well as other aspects 
of this GQTM are given elsewhere [14]. 

For this GQTM motion from state \m, i) to state \m', i) is potential free unless 
the 5th terms of T or T* are active. In this case, for \m,i) = \l m ,j m , S m {j m )) 
l m — 2 and S m (j m ) = 1 or(and) l m _ 1 = 2 and S'rn-iO'm-i) — 1- If either one 
but not both of these conditions hold, then by Eq, 5, V — K(l — 7). If both 
conditions hold then V = 2K(1 — 7). Equivalent conditions can be given for 
the state \m',i). 

An analysis of the types of potential barriers that can occur in a distribution 
for this shows that they are at least 2 path states wide. A barrier that is m + 1 
sites wide, corresponding to m successive (read 1) steps in which the 5th term 
of T or T' is active, has a central core of height V = 2K(1 — 7) occupying 
to — 1 successive path sites or states. The core is flanked by two potentials of 
height V = K(l — 7) occupying just one site or state. For to = 1 the barrier 
occupies 2 sites and is of height K(l — 7). 

For potential free regions the relation between energy and momentum is given 
by E = 2K(1 — cos k). For the core region of a barrier E = 2k(l — 7 cos /) (or 
E = 2K(1 — 7 cosh/) if E < V). From these equations one obtains cos A; = 
7 cos /(or cos k = 7 cosh /). For the flank regions E = K(2 — ( r y 2 + 27 cos 2h + 
l) 1 / 2 but this relation turns out not to be needed. 

The above describes all possible types of barriers that can occur on any com- 
putation path. However the distribution of barriers according to their width 
and spacing between barriers can vary widely among different paths. 

To apply this to a specific example, consider the initial state, shown in Figure 
1, with the head in state |0) in a wave packet localized to the left of the origin. 
All qubits are in state |0) except those at sites 0, A^ + 1 which are in state |2). 
The initial head and lattice qubit state are the path labels for this example. 
Iteration of T on this state generates in turn all the integers as binary strings 
of length < N . When the space between the two markers is completely filled 
with Is, corresponding to the integer 2^ — 1, the last pass of the head changes 
all Is to 0s. The head in state |1) then moves to the right away from the 
marker region as the enumeration is completed. A more detailed description 
of this process, based on iteration of T, is given elsewhere [13]. The potential 



Fig. 1. Initial and Final States for Counting GQTM for the First 2 n Binary Numbers. 
All lattice qubits are in state |0) except those at sites and n + 1 which are in state 
1 2). The initial and final head states are shown as wave packets with internal head 
states |0) and |1) to the left and right respectively. 

distribution associated with this example is obtained from a function IZn from 
the integers to {0, 1} such that TZisr{m) = 1 if and only if the 5th term of T, 
Eq. 10 is active in the transition T\m,i) — > \m + l,i) (i.e a read-1 step). The 
potential distribution, obtained from TZn and Eq. 5 [14] is given by 

(m'\Vi\m) = K(\ - 7)[Jl N (m)6 m i >m+1 + K N (m - l)5 m ', m -i]- (11) 

Note that for each component in the initial head wave packet state of Figure 
1 there is a different function TZn- This difference is ignored here because 
the functions differ only by a translation and describe the same potential 
distribution [13]. 

It turns out that the distributions of Is and 0s in TZn, exclusive of the 
tails, can be represented as an initial segment of length 2 N of the heirarchical 
sequence [16,17]. If one defines the integers m as finite strings of 0s and Is by 
m = 01 m m+1 where l m and m+1 denote strings of m Is and m + 1 0s, then 
(except for the tails) IZn = Q, 1, 0, 2, 0, 1, 0, 3, • • • , N_- This definition of TZn 
is based on the initial state [1, N + 1, S) with \S) as in Figure 1 chosen to be 
the path origin. 

This sequence is an initial segment of length 2 N of the heirarchical sequence 
in the underlined numbers. Based on this it can be shown that IZn can be 
generated by a pair of recursion relations given by [13,14] 



for n= 1,2, ...,JV with S = Q. 

Since the recursion relations generate the potential distribution it can be seen 
that the transfer matrix Zj^ is given by a pair of recursion relations obtained 
from Eq.12 as 

Z n = W n X n -i 

X n =X n _iW n X n -i (13) 



with Xo|n = e 2lk and X0I12 = for n — 1, 2 • • • , N. W n is the transfer matrix 
for the potential corresponding to the sequence n = 01™0™ +1 . Note that the 
order of matrix multiplication is the inverse of the order in which terms appear 
in the sequence of Eq. 12. 

The main advantage of these recursion relations is that only polynomially 
many matrix multiplications are needed to obtain the elements of the matrix 
Zjy. If Zjy is obtained from the matrices associated with each potential bar- 
rier in the distribution, then exponentially many matrix multiplications are 
required. 

From Eq. 13 one sees that the only matrices needed in explicit form (other 
than Xq) are the W m for m < N. These are given explicitly by [14] 



ik(m+2) 

W m \ u = [e 2lk sin lm - 2 7 e lk sin l(m - 1) + 7 2 sin l(m - 2)1 

2i 7 sin k sin 1 

pikm 

Wmli2 = -[sin/m — 27 sin /(m — 1) cos k + r ) 2 sin/(ra — 2)1.(14) 

2i7 sin k sin 1 

Here the momenta /, k, which refer to the core and potential free regions, are 
related by cos A; = 7 cos/. Note that W m \ n = W^\ 2 2, W m \ 2 i = W^\ 12 and W m 
is unimodular. 

These equations are valid for the unbound region E > V. The equivalent 
equations for the bound region are obtained by replacing sin everywhere in 
Eqs. 14 by isinh, the hyperbolic sine function. Also cos A; = 7 cosh/. 



5 Landauer Resistance, Band Spectra 



Calculations are made here for the LR and band spectra for several different 
values of 7 and iV and momentum regions for the counting GQTM. The band 
spectra refer to the associated periodic system which is a crystal whose unit 
cell is the potential distribution corresponding to enumeration of the first 
2^ integers. The momentum regions include the bound and unbound state 
regions. There are 7 dependent lower bounds on the values of \k\ below which 
there are no bands and the LR is very large [14]; these regions of k values are 
excluded from the calculations. 

Each figure gives a plot of the Log (base 10) of LR as a function of momen- 
tum k for chosen values of N and 7 and a region of k values. The log LR is 
plotted instead of the LR because the LR fluctuates rapidly and violently over 
many orders of magnitude and the interest here is in regions of appreciable 
transmission (LR values of order unity). The band and gap spectra for the 
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Fig. 2. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 10 and 7 = 0.999 for a Momentum Range 0.0223 < k < 0.10. 
The energy band spectrum for the associated periodic system is shown at the bottom 
of the figure where upper horizontal line segments correspond to energy bands and 
lower horizontal segments correspond to energy gaps. Very short bands or gaps 
appear as points. Band-gap edges show as vertical lines. The ordinal placement of 
the band spectrum, at the bottom of the figure, is done for convenience only. 

same parameters is also included so that one can compare the LR fluctuations 
to the band spectra. The number of momentum values used to generate the 
curves in each figure ranges from about 3,000 to 8,000. 

As both the LR and band spectra are extremely sensitive to values of 7, values 
that are reasonable physically should be chosen. This is done [14] by recalling 
that if T is the bilateral shift on the qubit lattice, the Feynman Hamiltonian. 
Eq. 1 is the symmetrized lattice equivalent of (— h /2m)d 2 /dx 2 . From this and 

V = 2Kb where 5 = 1—7, one obtains 5 = VmA 2 /h where A is the qubit 
lattice spacing. For electron systems, m is of the order of the electron mass, 

V is a few electron volts, and A is measured in Angstroms. Taking m equal 
to 2 electron masses, A = 1A, and V = 2ev gives S ~ 0.001 or 7 ~ 0.999. 

Using this value of 7 the log LR and band spectra for the associated periodic 
system have been calculated for N = 10 for several momentum or energy 
regions. The initial qubit state is shown in Figure 1. The results are shown in 
Figures 2 and 3 for momentum ranges of 0.0223 < k < 0.10 and 0.095 < k < 
0.220. The bound state region extends up to k = 0.0447 with the unbound 
region extending up and beyond the upper limit. In this and all succeeding 
figures the band spectrum is shown as a band with upper or lower flat line 
segments denoting bands or gaps respectively. Very narrow gaps (or bands) 
are shown as downward (or upward) pointing spikes. The results in the two 
figures show that the characteristics of the band spectra and log LR found 
for k < 0.0447 [14], extend smoothly into the region where E > V. Most 
of the region is occupied by bands with most gaps being quite narrow. A 
few wider gaps, such as those at k — 0.033, 0.077, 0.149 and especially at 
k = 0.054, 0.101, 0.198 are present. 

The log LR can be characterized as a downward trending band of fluctuating 
values from an average value of about — 1, (LR = 0.1) at the lower end of the 
k region to about —2, (LR = 0.01) at the upper end of Figure 2 to a value of 
about —3.5 at the upper end of Figure 3. Several peaks extending to higher 
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Fig. 3. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 10 and 7 = 0.999 for a Momentum Range 0.095 < k < 0.220. 
The energy band spectrum for the associated periodic system is shown at the bottom 
of the figure. Additional details are given in the caption for Figure 2. 



Fig. 4. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 18 and 7 = 0.999 for a Momentum Range 0.0223 < k < 0.10. 
The energy band spectrum for the associated periodic system is shown at the bottom 
of the figure. Additional details are given in the caption for Figure 2. 

values (up to LR ~ 10) project out of the band. 

The peaks in log LR appear to be associated with gaps, with the higher peaks 
associated with wider gaps. This effect [14] was also found for calculations 
with the Kronig- Penney model with potential distributions corresponding to 
the Fibonacci and Thue- Morse substitution sequences [20,21]. The minima in 
log LR which show as downward pointing spikes appear to be associated also 
with narrow gaps. Additional details on this are given elsewhere [14]. 

Results of calculations of the band spectrum and log LR for a larger value of 
A" = 18, corresponding to a potential distribution with 2 17 individual potential 
barriers, are shown in the next two figures for the same two momentum ranges 
that were used in Figures 2 and 3. The same value of 7 = 0.999 was used. The 
band spectra are much more finely divided with each relatively broad band for 
N = 10 divided into a great number of very narrow bands with intervening 
narrow gaps. The fraction of the momentum regions occupied by bands is 
smaller than that for A" = 10. The LR fluctuates very rapidly over many orders 
of magnitude between low values less than 10 -4 to high values greater than 
10 20 . For much of the momentum region log LR occupies a band (the black 
region in the figures) which trends down from a value around (LR = 1) at 
the low momentum end of Figure 4 to a value around —3 (LR = 10 -3 ) at 
the high momentum end of Figure 5. The trend, which represents a kind of 
average, is in the right direction as the LR should decrease as the energy is 
increased. 
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Fig. 5. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 18 and 7 = 0.999 for a Momentum Range 0.095 < k < 0.22. 
The energy band spectrum for the associated periodic system is shown at the bottom 
of the figure. Additional details are given in the caption for Figure 2. 



Fig. 6. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 10 and 7 = 0.99 for a Momentum Range 0.069 < k < 0.22. The 
energy band spectrum for the associated periodic system is shown at the bottom of 
the figure. Additional details are given in the caption for Figure 2. 

Probably the most remarkable aspect is that the transmission probability is 
high (> 50 — 99%) on average, yet it is filled with numerous holes as small 
regions of very low or negligible transmission. This is shown by the numerous 
excursions of log LR to high values > 4 from low average values around to 
—2 and lower in Figures 4 and 5. It is noteworthy that this effect appears to be 
independent of whether E<VorE>V(E = Vatk = arccos7 = 0.0447). 
Also there is a correlation between the widths of regions of low transmission 
and the associated band spectrum gap widths in that wider gaps are correlated 
with wider regions of low transmission. Examples are the relatively wide band 
gaps at k = 0.054, 0.100, 0.198. For these gaps the widths of the regions of 
very low transmission (high values of log LR) are relatively wide. 

Calculations made for the same values of A^ = 18 and 7 = 0.999 for higher 
values of the momentum, up to k = 0.5, show that the pattern described 
above continues. The "black" band containing much of the fluctuating values 
of log LR continues a slow downward trend to f» —4 at the upper end of the 
region. Numerous narow spikes to high values are present. In some narrow 
momentum regions the spikes are so dense that high resolution calculations 
are needed for their individual resolution. 

In order to investigate the dependence of the LR and band spectra on 7, 
calculations were made for a value of S larger by a factor of 10, i.e. 7 = 0.99, 
and A" = 10. The results for the energy regions 0.069 < k < 0.220 and 
0.20 < k < 0.42 are shown in the next two figures. 
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Fig. 7. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 10 and 7 = 0.99 for a Momentum Range 0.20 < k < 0.42. The 
energy band spectrum for the associated periodic system is shown at the bottom of 
the figure. Additional details are given in the caption for Figure 2. 

The band spectra show the presence of a number of wide gaps. Also most 
of the bands are quite narrow. Comparison with the band spectra of Figs. 2 
and 3, which are for the same value of N but 5 smaller by a factor of 10, 
show that the fraction of the momentum regions occupied by bands is much 
lower and gap widths are wider, log LR shows a band of fluctuating values 
decreasing slowly from an average value of ~ 1(LR = 10) at the low end of 
the momentum region in Figure 6 to a value of ~ —2(LR = 0.01) at the upper 
end of the momentum region in Figure 7. 

Superimposed on this are peaks in the values of log LR which extend up to 
18 or more. Each peak appears to occupy the same momentum region as the 
gap, with the peak height correlated with the gap width in that wider gaps 
are associated with higher peaks. This is especially clear in Figure 6. The 
distribution of peak heights in the figure also shows a heirarchical relationship 
similar to that shown by the potential distribution of Eqs. 11 and 12. That 
is, the region between the two largest peaks is divided in two by a peak of 
lesser height. Each subregion is divided in two by an even smaller peak. This 
division into two subregions by peaks of decreasing height continues down to 
the smallest peaks associated with the small gaps and which form the band 
of fluctuating values. This regularity is less evident in the momentum region 
shown in Figure 7 for which the fluctuations in log LR are more chaotic. Since 
the value of k = 0.142 for which E — V occurs in the middle of Figure 6, 
it appears that the regularity of the peak distribution is not dependent on 
whether the momentum is in the bound or unbound region. 

In order to investigate in more detail the relation between peaks in log LR and 
bands and gaps, a high resolution calculation for 7 = 0.99, A/" = 10 was done 
for the small momentum region 0.21 < k < 0.23. The results are shown in 
the next figure. The results show that each peak in log LR is defined by sharp 
minima which are associated with relatively narrow gaps. More exactly the 
minima seem to be associated with the band-above-gap edges of the narrower 
gaps. In this fashion the momentum region of each peak includes one band and 
one gap. Each gap has a sharp minima associated with its upper edge except 
for the wider gaps. This appears to result from the fact that in general the 
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Fig. 8. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 10 and 7 = 0.99 for a Small Momentum Range 0.21 < k < 0.23. 
The energy band spectrum for the associated periodic system is shown at the bottom 
of the figure. Additional details are given in the caption for Figure 2. 



Fig. 9. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 8 and 7 = 0.9 for a Momentum Range 0.18 < k < 0.52. The 
energy band spectrum for the associated periodic system is shown at the bottom of 
the figure. Additional details are given in the caption for Figure 2. 



Fig. 10. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 8 and 7 = 0.9 for a Momentum Range 0.5 < k < 1.0. The 
energy band spectrum for the associated periodic system is shown at the bottom of 
the figure. Additional details are given in the caption for Figure 2. 

depth of the minima in log LR decreases with increasing width of the gap at 
whose edge it is located until it disappears entirely for sufficiently wide gaps. 
For these wider gaps the momentum region of the associated peak includes 
two bands and two gaps. This effect, which was noted before [14], is seen for 
the wide gaps centered at k = 0.216, 0.2185, 0.221, 0.2267. 

Increasing S by another factor of 10 to 0.1 with 7 = 0.9 gives the results 
shown in the next two figures for the momentum regions 0.18 < k < 0.52 and 
0.5 < k < 1.0 respectively. The calculations are for N = 8. Since k = 0.45 for 
E — V, the energy region used in Figure 9 is mostly the bound state region. 
That used in Figure 10 includes part of the unbound state region. The results 
in the figures show that the potential is high enough (10 times the height 
of that in Figures 6 and 7) so that definite differences in the band spectra 
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and log LR are evident for the two regions. In the bound state region (Fig. 
9) the bands are extremely narrow and are few and widely separated with 
many large gaps present. The values of log LR show the the LR is very high 
for most of the region, except for a few very narrow spike minima to values 
around 1, (log Li? = 0). The spike minima occur at the same energies as the 
bands. 

In the unbound region (Fig. 10) the bands are much more numerous and in- 
dividual bands are wider than in the bound state region. The LR is much 
lower on average, and there are extended regions (0.59 < k < 0.62 and 
0.70 < k < 0.81) for which LR ~ 1. To summarize, the results for these 
figures show that in the bound state region there is almost no transmission 
through the 128 potentials except for a few extremely narrow regions at which 
bands are located. In the unbound state region there are wide regions of good 
transmission of ~ 50% or more. Energy bands are more numerous and are 
wider. 



6 discussion 



The values of 7 and N for which results of calculations of band spectra and 
log LR have been shown are quite limited. These limitations are imposed by 
the requirement that the calcuations of the band spectra and LR be reasonably 
reliable. This is not trivial because fluctuations in the values of the matrix ele- 
ments of Z^ used to calcuate the LR and the band spectra become extremely 
rapid and extreme at higher values of iV and lower values of 7. For example, 
even for N = 18 and 7 = 0.999 or N = 8 and 7 = 0.9, TrZ^ and the LR 
fluctuate rapidly between ±10 M with M taking values of 200-300 or more for 
some values of k. Also these fluctuations occur over extremely small intervals 
of A;. 

In essence limitations on the values of N and 7 are imposed by (classical) 
computer program limitations on the magnitude and number of significant 
figures allowed in the computations. Modifications of the programs to allow 
calculations over a larger range of iV and 7 values are not warranted as they 
are complex and at most give one a slight extension of the range of acceptable 
values of N and 7. 

In many ways the results obtained here, which are similar to those obtained by 
Roy and Khan [20,21], are unexpected. The LR fluctuates rapidly with many 
spikes down to extremely low values over the momentum range examined. 
Except for the highest potential considered (i.e. 7 = 0.9) there is not much 
difference between the bound and unbound regions as far as the band spectra 
or log LR are concerned. The transmission, even in the bound region for N = 
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18 and 7 = 0.999, is high (~ 50% or more) on average although there are a 
great many holes of extremely low transmission. 

These results need explaining since one would expect little or no transmission 
in the bound region, because of decay in the transmission amplitude resulting 
from tunnelling effects and localization [10]. This is especially the case for 
N = 18 with more than 128, 000 potential barriers present. Also one would 
expect almost complete transmission in the unbound regions with few or no 
regions of low transmission. Neither of these results occur. 

It is suggested here that the calculated results are best understood in terms of 
coherent superposition of the many reflected waves moving backwards through 
the potential distribution. At some values of the momenta the reflected waves 
interfere constructively, which greatly increases the amplitude of the overall re- 
flected component and decreases the amplitude to the transmitted component. 
At other momentum values the reflected waves interfere destructively. In this 
case the amplitude of the reflected component is small and the amplitude of 
the transmitted component is high. Q Note that if B and F are the complex 
valued reflection and transmission coefficients respectively for the potential 
distribution, Eqs. 6 and 7, then normalization requires that |-B| 2 + |-F| 2 = 1. 

The number of component waves which add coherently should be very large. 
Since each barrier in the distribution generates transmitted and reflected com- 
ponents from any impinging wave, whether it is the input wave, or waves re- 
flected or transmitted from other barriers, the number of multiple reflections 
or transmissions or combinations of the two becomes extremely large, espe- 
cially in a distribution with many barriers. All the components moving in the 
same direction add coherently. 

This coherent superposition is expected to depend sensitively on all parame- 
ters involved. It also should be very dependent on the fact that the potential 
distribution is not random; in the example considered here the distribution is 
quasiperiodic [24,25] and is heirarchical. Intuitively if coherent superposition 
is the dominant effect, one would expect that for fixed potential height (fixed 
7) , the fluctuation rate of the Landauer Resistance and the number of bands 
and gaps to increase with increasing N. The reason is that the number of in- 
dividual potential barriers present increases exponentially with increasing N. 
The corresponding increase in the number of reflected and transmitted compo- 
nent waves which combine coherently means that momentum regions in which 
backward moving reflections interfere destructively to give high transmission 



1 Exactly the same argument applies mutatis mutandis. Constructive or destructive 
interference of the transmitted components decreases or increases the amplitude of 
the reflected component emerging from the distribution. From the viewpoint of a 
computation one is more interested in the amplitude of the transmitted component 
than in the reflected one. 
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Fig. 11. The magnitudes of the Reflection (B) and Transmission (F) Coefficients for 
a Single Potential Barrier Corresponding to m = 10 read 1 Steps for 7 = 0.999 for 
a Momentum Range of < k < 0.5. The values of \B\ for k > 0.3 were increased by 
a factor of 10 to show more clearly the nodal structure. 

become much narrower and broad bands split into many narrower bands with 
intervening gaps. 

In a similar fashion momentum regions in which backward moving reflections 
combine constructively to give very low transmission also become narrower. 
This narrowing and splitting of regions with constructive or destructive inter- 
ference results in a large increase in the extent and rapidity of fluctuations in 
the LR and splitting of the energy bands. 

The importance of coherent effects in the transmission and reflection of a state 
through the potential distribution is demonstrated clearly by consideration of 
the transmission and reflection for a single potential barrier. For motion in 
the potential distribution corresponding to iteration of all binary numbers of 
length < N the complex transmission (F) and reflection (B) coefficients are 
obtained from Eqs. 7 and 8 and the normalization |-F| 2 + \B\ 2 = 1 as 

F = 77V, B = -I4 1 - (15) 



J N 22 ^N 22 



For a single barrier corresponding to m read 1 steps these equations hold with 
Zn replaced by W m , Eq. 14. The momentum dependence of the magnitudes 
and phases of F and B for m = 10 and 7 = 0.999 for a momentum range of 
< k < 0.5 are shown in the next two figures. The magnitude values for the 
reflection coefficient B for k > 0.3 have been multiplied by a factor of 10 to 
illustrate the nodal structure of B more clearly. The range of the phases in 
Figure 12 extends from — it to n. The large discontinuities in the curves for 
both B and F are a artifact of the presentation method in that tt + x = — n + x. 
The discontinuity in the phase for B at k ~ 0.3 is real, though. It occurs at 
the node point for \B\, Figure 11. The figures show that the magnitudes of 
B and F tend rapidly to values close to and 1. At a value of k — 0.0447 
corresponding to E + V, \F\ = 0.976 and \B\ = 0.218. Thus at momenta 
corresponding to energies above the barrier height there is still appreciable 
reflection. Unlike the magnitudes the phases of both B and F vary steadily 
over the range of momenta shown. This is significant because it means that 



Fig. 12. The Phases of the Reflection (B) and Transmission (F) Coefficients for a 
Single Potential Barrier Corresponding to m = 10 read 1 Steps for 7 = 0.999 for a 
Phase Range of [— ir, n]. The large discontinuities in the phases for B and F are a 
result of the figure presentation in which —ir and ir are not identified. 

coherent additions of the many reflected and transmitted components in a 
sequence of these barriers would be expected to give rapidly varying overall 
transmission and reflection coefficients. This is shown by the calculations. 

Calculations show that as N is increased for constant 7 the internodal maxima 
increase and the distances between nodes in \B\ decrease. The rate of change 
of the phases of both F and B also increases and the discontinuities at the B 
node points increase. This implies an increase in the rapidity of fluctuations 
of the LR and in the number of spectral bands and gaps with increasing N, 
also in agreement with the calculations. 

As the potential is increased (7 decreased) but iV is fixed, the position of the 
first node in \B\ moves out to higher values of k. Also the internodal maxima 
increase. However the rates of change of the phases for both B and F decrease 
out to the first node as the potential height increases. For higher values of 
k the phase slopes appear to be roughly independent of the potential height. 
It is of interest to note that the position of the first node at k — 0.59 for 
7 = 0.9, N = 8 is close to the value of k at which the LR and band spectra 
change character (See discussion of Figures 9 and 10). 

The curves of Figure 11 can be used to show the importance of coherent effects. 
Consider for example the potential distribution for N = 18 and 7 = 0.999 used 
to calculate the LR and band spectra shown in Figures 4 and 5. The proper- 
ties of the heirarchical distribution are such that for m — 1, 2, • • • , N — 1 there 
are 2 N ~ m ~ 1 potential barriers of width corresponding to m read 1 steps in the 
distribution. For the example at hand there are 128 barriers for m = 10. If the 
reflections and transmissions from these barriers combine incoherently, then 
the probability of transmission through just these barriers and ignoring the 
others in the distribution is given by (1 — |5(A;)| 2 ) 128 where the k dependence 
is shown explicitly. It is clear from the k dependence of |-B(fc)l shown in Figure 
12 and the LR shown in Figures 4 and 5 that the k dependence of the prob- 
ability of transmission is completely different from the incoherent prediction. 
Inclusion of the effects of the other barriers for different m in the incoherent 
prediction whould make the disagreement even more extreme. 
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The coherent effects implied in the results for the LR and band spectra give 
a good illustration of nonlocality in quantum mechanics. Consider a value 
of k at which the reflection components interfere destructively as they move 
backwards along the path. This means the amplitude of the overall reflection 
coefficient for the reflected wave emerging from the distribution is a (local) 
minimum and the amplitude for the overall emerging transmitted wave is a 
maximum. This is a nonlocal effect since much of what happens to the reflected 
components occurs at path locations distant from those of the transmitted 
components. The effect is quite pronounced since it extends over the whole 
potential distribution. In the case of N — 18 this effect extends over 10 6 path 
sites. The distance over which this effect occurs grows exponentially with N as 
the number of path sites in the potential distribution is about equal to 2 N+2 . 

Is is to be emphasized that this effect is quite general and is not restricted 
to the specific distribution considered here. It is shown, for instance, in the 
graphs of the LR versus energy for different potential distributions obtained 
by Roy and Kahn [20,21]. These effects also play an important role in most 
work on transmission of electrons through nonrandom potential distributions. 
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